Studying quantum spin systems through entanglement estimators 
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We study the field dependence of the entanglement of formation in anisotropic 5=1/2 antiferro- 
magnetic chains displaying a T — field-driven quantum phase transition. The analysis is carried 
out via Quantum Monte Carlo simulations. At zero temperature the entanglement estimators show 
abrupt changes at and around criticality, vanishing below the critical field, in correspondence with 
an exactly factorized state, and then immediately recovering a finite value upon passing through 
the quantum phase transition. At the quantum critical point, a deep minimum in the pairwise-to- 
global entanglement ratio shows that multi-spin entanglement is strongly enhanced; moreover this 
signature represents a novel way of detecting the quantum phase transition of the system, relying 
entirely on entanglement estimators. 
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PACS numbers: 03.67.Mn, 75.10.Jm, 73.43.Nq, 05.30.-d 

Collective behavior in many-body quantum systems 
is associated with the developemcnt of classical corre- 
lations, as well as of correlations which cannot be ac- 
counted for in terms of classical physics, namely entangle- 
ment. Entanglement represents in essence the impossibil- 
ity of giving a local description of a many-body quantum 
state. In particular entanglement is expected to play an 
essential role at quantum phase transitions, where quan- 
tum fluctuations manifest themselves at all length scales. 
The behavior of entanglement at quantum phase transi- 
tions is a very recent topiCjSO far investigated in a few ex- 
actly solvable cases 0,12 00 Moreover, entanglement 
overwhelmingly comes into play in quantum computation 
and communication theory, being the main physical re- 
source needed for their specific tasks ||. In this respect, 
the perspective of manipulating entanglement by tunable 
quantum many-body effects appears very intriguing. 

In this letter we show that entaglement estimators give 
important insight in the physics of spin systems. In par- 
ticular, we focus on two striking features of anisotropic 
spin chains in an external field: the occurrence of a fac- 
torized ground state at a field hf and of a quantum phase 
transition at h c . We propose a novel estimator to under- 
stand the role of quantum fluctuations in the quantum 
critical region. 

We focus our attention on the ID XYZ model in a field: 
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where J > is the exchange coupling, i runs over the sites 
of the chain, and h = g^H/J is the reduced magnetic 
field. In Eq. Q we have implicitly performed the canon- 
ical transformation S*' y — ► (— iySf' v with respect to the 
more standard antiferromagnetic hamiltonian. The pa- 
rameters A y , A z > control the anisotropy of the sys- 
tem. In particular, for A z = Eq. reduces to the 



exactly solvable XY model in a transverse field 0. For 
A z 7^ the model does not admit an exact solution 0, 
and it has been recently studied within approximate ana- 
lytical and numerical approaches H,E1E3- Interestingly, 
the general model with finite A z is experimentally real- 
ized by the 5 = 1/2 quantum spin chain CS2C0CI4 [llj . 
with strong planar XZ anisotropy, A y rs 0.25, A z 1, 
and J sa 0.23 meV. 

In our study, we concentrate on the case < A y < 1, 
A z = 1, defining the XYX model in a field 12] . The more 
general case of the XYZ model in a magnetic field with 
A z < 1 should exhibit the same qualitative behavior, 
as it shares the same symmetries with the XYX model 
in a field, and therefore it belongs to the same univer- 
sality class. The analysis is carried out via Stochastic 
Series Expansion (SSE) Quantum Monte Carlo (QMC) 
simulations, based on a modified directed-loop algorithm 
[T3L fl4| . for chains of various length, from L = 40 to 
L = 120. Ground state properties have been determined 
by considering inverse temperatures /3 = 2L, in order to 
capture the T = behavior. 

The ground-state phase diagram of the XYX model in 
the A y — h plane is shown in Fig.^ The model displays 
a field-driven quantum phase transition at a critical field 
h c (A y ), which separates the Neel-ordered phase (h < 
h c ) from a disordered phase (h > h r ) with short-range 
antiferromagnetic correlations 0, El EjJ • The transition 
line h c (A y ) has been determined by a scaling analysis 
of the correlation length £ xx , whose linear scaling £ xx ~ 
L marks the quantum critical point. Using the critical 
scaling of the structure factor S xx (q = 0) ~ L 1 / U ^ z , we 
verified that the transition belongs to the ID transverse- 
field Ising universality class (7/f = 7/4 and z — 1), in 
agreement with analytical predictions 

Besides its quantum critical behavior, a striking fea- 
ture of the model of Eq. JJJ is the occurrence of an 
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FIG. 1: Upper panel: ground state phase diagram of the 
XYX model in a field. Mean-field (MFA) results are taken 
from Ref. Iinl the factorizing field h{ from Ref. 0. Lower 
panel: quantum critical behavior of a;— and z— magnetizations 
and correlation length for the model with A y — 0.25, L — 100, 
P — 200. The factorizing field is indicated by a dashed line. 
The arrow indicates the quantum critical point (QCP). 



exactly factorized ground state for a field hf (A y Mower 
than the critical field h c , as predicted in Ref. [jj. In 
the case of the XYX model, this factorizing field is Q 
hf = y2(l + A y ). At h = hf the ground state of the 
model takes a product form |W) = ® i=1 where the 
single-spin states \ipi) are eigenstates of («i(2) • S) with 
n i(2) being the local spin orientation on sublatticc 1 (2). 
Taking n = (cos</>sin(9,sin^>sin#,cos(9), one obtains [j| 

0, 4> 2 = 7T, 



, sin sin a, co 

cos" 1 



one obtains 

The fac- 



>l = o, cp 2 — n, ui = (72 = cos y^i-r ^»)/2 
torized state of the anisotropic model continuously con- 
nects with the fully polarized state of the isotropic model 
in a field for A y — 1 and h = 2. 

Despite its exceptional character, the occurrence of a 
factorized state is not marked by any particular anomaly 
in the experimentally measurable thermodynamic quan- 
tities shown in the lower panel of Fig. ^ However, we 
will now see how entanglement estimators are able to pin 
down a factorized state with high accuracy. 

To estimate the entanglement of formation [l5( in 
the quantum spin chain of Eq. JTJ) we make use of 
the one-tangle and of the concurrence. The one-tangle 
[HI E3 quantifies the T = entanglement of a sin- 
gle spin with the rest of the system. It is defined as 
n = AdetpW, where = (I + J2 a M a a a )/2 is the 
one-site reduced density matrix, M a = (S a ) (estimated 
as M a = |(Sf% L/2 )| 1/2 in the QMC on a chain of 
length L), o~ a are the Pauli matrices, and a — x,y,z. 
In terms of the spin expectation values M a , n takes the 
simple form: 

r 1 = l-4]T(M*) 2 . (2) 

a 

It can be easily shown that the vanishing of t\ implies a 
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FIG. 2: (a) One-tangle ri and sum of squared concurrences ti 
as a function of the applied field for the S = 1/2 XYX model 
with A y = 0.25, L = 100 and /3 = 200. Inset: contributions to 
the concurrence between j-th neighbors. Full symbols stand 
for Cj-V+j, open symbols for C^ +J -. The dashed line marks 
the critical field h c . (b) Concurrence as a function of spin-spin 
distance. Parameters as in the previous figure. 



factorized ground state, and viceversa. 
The concurrence [lij quantifies instead the pairwise en- 
tanglement between two spins at sites i, j both at zero 
and finite temperature. For the model of interest, in ab- 
sence of spontaneous symmetry breaking (M x — 0) the 
concurrence takes the form |l7j 



C l3 -2max{0,Cy ) ,G^ ) } , (3) 



where 



C$ = 



ci 2) = 



\9??+9l?\ 




(4) 



(M*)2, (5) 



with g%* = {SfSf). 

The T = QMC results for the model Eq. JT} with 
A y = 0.25 are shown in Fig. |2 where we plot ri, the 



sum of squared concurrences 



T2 



(6) 



and the field and space dependence of the concurrence. 
The following discussion, although directly referred to 
the results for A y = 0.25, is actually quite general and 
applies to all the other studied values of A y . 
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FIG. 3: Entanglement ratio t^/ti as a function of the field 
for Aj, = 0.25 and /3 = 2L. Inset: entanglement ratio for the 
ID Ising model in a transverse field. The dashed line is the 
L — > oo extrapolation result. 



Unlike the standard magnetic observables plotted in 
Fig. ^ the entanglement estimators display a marked 
anomaly at the factorizing field, where they clearly van- 
ish as expected for a factorized state. When the field is 
increased above /if, the ground-state entanglement has a 
steep recovery, accompanied by the quantum phase tran- 
sition at h c > h{. For A y — 0.25, e.g., h c = 1.605(3) and 
h{ = 1.5811. The system realizes therefore an interest- 
ing entanglement switch effect controlled by the magnetic 
field. 

As for the concurrence, Fig. |2Ib) shows that its range 
is always finite at and around the critical point, and it 
never extends farther than the fourth neighbor. More- 
over, the factorizing field divides two field regions with 
different expressions for the concurrence: 

C\f < < C\f for h<ht, (7) 
C\f < < C\f for h > hi , (8) 

whereas c[) ] = df = at h = h{. 

In presence of spontaneous symmetry breaking occur- 
ring for h < h c , the expression of the concurrence is gen- 
erally expected to change with respect to Eqs. Q,©, 
as extensively discussed in Rcf. For the model un- 

der investigation, this happens when the condition [l9| 

df < C$ is satisfied, i.e. for h > h { . This means 
that our estimated concurrence is accurate even in the 
ordered phase above the factorizing field; in the region 
< h < hi it represents instead a lower bound to the 
actual T = concurrence. Alternatively it can be re- 
garded as the concurrence for infinitesimally small but 
finite temperature. 

In Fig. |2 the sum of squared concurrences T2 is always 
smaller than or equal to the one-tangle t\ , in agreement 
with the Coffman-Kundu-Wootters conjecture 16] . This 
shows that entanglement is only partially stored in two- 
spin correlations, and it is present also at the level of 
three-spin entanglement, four-spin entanglement, etc. In 
particular, we interpret the entanglement ratio R = t%/ti 
as a measure of the fraction of the total entanglement 



stored in pairwise correlations. This ratio is plotted as a 
function of the field in Fig. [3J As the field increases, we 
observe the general trend of pairwise entanglement sat- 
urating the whole entanglement content of the system. 
But a striking anomaly occurs at the quantum critical 
field h c , where R displays a very narrow dip. According 
to our interpretation, this result shows that the weight 
of pairwise entanglement decreases dramatically at the 
quantum critical point in favour of multi-spin entangle- 
ment. Unlike classical correlations, entanglement shows 
the special property of monogamy 16] , namely full entan- 
glement between two partners implies the absence of en- 
tanglement with the rest of the system. Therefore multi- 
spin entanglement appears as the only possible quantum 
counterpart to long-range spin-spin correlations occur- 
ring at a quantum phase transition. This also explains 
the somewhat puzzling result that the concurrence re- 
mains short-ranged at a quantum phase transition (Fig. 
Hfb)) while the spin-spin correlators become long-ranged, 
and it evidences the serious limitations of concurrence 
as an estimate of entanglement at a quantum critical 
point. Strong indications on the relevance of multi-spin 
entanglement in quantum-critical spin chains come also 
from the study of the entanglement between a block of 
L contiguous spins and the rest of the chainQ- Finally, 
multi-spin entanglement involves a macroscopically co- 
herent superposition of quantum states, and this result 
is consistent with the picture of macroscopic (i.e. long- 
wavelength) quantum fluctuations occurring at a quan- 
tum phase transition. 

In turn, we propose the minimum of the entanglement 
ratio R as a novel estimator of the quantum critical point, 
fully based on entanglement quantifiers. This result ap- 
pears general for the whole class of models described by 
the hamiltonian of Eq. Q . Inset (b) of Fig. [3] shows in 
fact that an analogous dip in the entanglement ratio sig- 
nals the quantum phase transition in the Ising model in a 
transverse field (A y = A z = 0), occurring at the critical 
field h c = 1/2. Work is in progress to test the universal- 
ity of such novel signature of quantum critical behavior 
for completely independent quantum phase transitions. 

The use of the QMC method enables us to naturally 
monitor the fate of entanglement when the temperature 
is raised above zero. In this regime the concurrence is the 
only well-defined estimator of entanglement, whereas the 
one-tangle has not yet received a finite-temperature gen- 
eralization. Fig. 01a) shows the nearest-neighbor (n.n.) 
concurrence as a function of the field for different temper- 
atures in the XYX model with A y — 0. We observe that 
Ci^+ n < C^i+i for n > 1, and at high enough tempera- 
ture (T > 0.1J) only the n.n. concurrence survives. The 
most prominent feature is the persistence of a field value 
(or an interval of values) at which the concurrence is ei- 
ther zero (for T > 0.05 J) or - 10~ 3 (for < T < 0.05 J). 
In particular, the field values for which the concurrence 
vanishes are temperature-dependent, so that two-spin en- 
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FIG. 4: (a) Nearest-neighbor concurrence Ci^+i as a func- 
tion of field for different temperatures in the XYX model 
with A y — 0, L = 40. One division on the y axis corre- 
sponds to a concurrence interval of 0.1. Offsets (marked by 
the dashed lines) have been introduced to improve readabil- 
ity, (b) Components of the nearest-neighbor concurrence and 
spin-spin correlators as a function of temperature at the fac- 
torizing field hf — v2- The dashed lines mark the zero value. 



mators, unlike the more conventional magnetic observ- 
ables. Moreover we find that entanglement estimators 
are able to detect the quantum critical point, marked by 
a narrow dip in the pairwise-to-global entanglement ra- 
tio. Therefore we have shown that entanglement estima- 
tors provide precious insight in the ground-state prop- 
erties of lattice S — 1/2 spin systems. Thanks to the 
versatility of QMC, the same approach can be used for 
higher-dimensional systems. In this respect, investiga- 
tions of the occurrence of factorized states in more than 
one dimension are currently in progress. Finally, the 
proximity of a quantum critical point to the factorized 
state of the system gives rise to an interesting field- 
driven entanglement-switch effect. This demonstrates 
that many-body effects driven by a macroscopic field are 
a powerful tool for the control of the microscopic entan- 
glement in a multi-qubit system, and stand as a profitable 
resource for quantum computing devices. 
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ing, G. Falci, R. Fazio, A. Osterloh, and G. Vidal are 
gratefully acknowledged. We acknowledge support by 
NSF under grant DMR-0089882 (T.R. and S.H.), by 
INFN, INFM, and MIUR-COFIN2002 (A.F., P.V., and 
V.T.). 



tanglemcnt can be switched on and off tuning both the 
field and the temperature. Fig. Hfb) shows a highly 
non-trivial temperature dependence of the two-spin en- 
tanglement at h = hf = y/2. Although vanishing at 
T = 0, the n.n. concurrence has a quick thermal acti- 
vation due to thermal mixing of the factorized ground 
state with entangled excited states. Although the spec- 
trum over the ground state displays a gap of order 0.1J 
H, in one dimension strong fluctuations induce thermal 
entanglement [2(1 ] already at temperatures which are an 
order of magnitude lower. The appearence of thermal en- 
tanglement is directly related to the increasing behavior 
of the correlators g vy = gff +1 and g zz = gf f +1 enter- 
ing the expression of the component, Eq. (also 
shown in Fig. EJb)). In particular the appearence of a 
finite g yy is a purely quantum effect, since A y = 0. Be- 
cause of the non- monotonous behavior of the correlators, 
at higher temperatures thermal entanglement disappears 
and reappears again, revealing an intermediate tempera- 
ture region where two-spin entanglement is absent. 

In summary, making use of efficient QMC techniques 
we have provided a comprehensive picture of the entan- 
glement properties in a class of anisotropic spin chains of 
relevance to experimental compounds. We have shown 
that the occurrence of a classical factorized state in these 
systems is remarkably singled out by entanglement esti- 
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